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Abstract 

Neutrino\antineutrino deep inelastic scattering (DIS) data provide useful constrains 
for the flavor decomposition in global flts of parton distribution functions (PDF). The 
smallness of the cross-sections requires the use of nuclear targets in the experimental 
setup. Understanding the nuclear corrections is, for this reason, of utmost importance 
for a precise determination of the PDFs. Here, we explore the nuclear effects in the 
neutrino\antineutrino-nucleon DIS by comparing the NuTeV, CDHSW, and CHORUS 
cross-sections to the predictions derived from the latest parton distribution functions 
and their nuclear modiflcations. We obtain a good description of these data and find no 
apparent disagreement between the nuclear effects in neutrino DIS and those in charged 
lepton DIS. These results also indicate that further improvements in the knowledge of 
the nuclear PDFs could be obtained by a more extensive use of these sets of neutrino 
data. 
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1 Introduction 



The undisputed success of the recent global analyses of the free nucleon parton dis- 
tribution functions (PDFs) [H 121 13] have proven the QCD factorization theorem as a 
well-working tool to analyse and interpret the experimental cross-sections measured 
at high-energy colliders. These analyses look for a universal set of PDFs, ff^{x,Q^), 
facilitating the cross-sections computations schematically as 

aAB^h+x = Yl ff(Q') ® ff(Q') ® ^ij-^h+x, (1) 

where &ij^h+x are perturbatively computable coefficients. 

The usual procedure to obtain a global fit of PDFs is well establishecl^: a functional 
form with parameters is assumed for the different PDFs at a given initial scale 
Qo- This initial condition is then evolved by the Dokshitzer-Gribov-Lipatov-Altarelli- 
Parisi (DGLAP) equations [1] and the best set of parameters obtained in an iterative 
procedure involving a quality criterion — a particular definition of the goodness of the 
fit through a generalized x^- Several differences among the analyses, like the choice of 
the initial conditions, the choice of the data sets, the actual definition of the or the 
approximations in the computation of the actual cross-sections, result in differences 
in the obtained PDFs. In particular, the neutrino\antineutrino DIS data is included 
only in part of these analyses due to the difficulty in removing the nuclear effects in 
a model-independent manner. These data are, however, of importance for the flavor 
decomposition of the PDFs which, on the other hand, affects the precision in the 
calculation of relevant cross-sections in present accelerators, in particular the LHC 

In the situation, where the nucleons are part of a bound nucleus, the factorization 
theorem could be theoretically more doubtful, as other processes like multiple scatter- 
ing could ruin its applicability. However, the corresponding global analyses O O [71 [8] 
show also an excellent description of charged-lepton deep inelastic scattering (DIS) 
and Drell-Yan (DY) dilepton production data involving nuclear targets with the as- 
sumption of universal, process-independent, nuclear PDFs (nPDFs). This universality 
has also been checked in other processes, as single-inclusive hadron production in dAu 
collisions at RHIC, flrst introduced in [9] and then also in [5]. The access to harder 
and harder probes in nuclear collisions made the need of this type of analyses, and the 
corresponding checks of Eq. ([T]), also of utmost importance for the phenomenological 
interpretation of the imminent LHC heavy-ion program. 

The amount of experimental information used in global nPDF flts is much smaller 
than the corresponding one for the free proton case. The total number of data points is 
(9(1000) scanning the range oi x > 10~^ in the perturbative region for different nuclei. 
These data are normally given in terms of the ratio of nuclear over free proton cross- 
sections (with the appropriate normalization factor) so most of the global analyses 
study the corresponding ratios of PDFs Rf{x,Q^) = /^°"'''^(x, g2)//frcc^^^ g2^)^ 

"^A new approach with a different treatment of the initial conditions, especially designed to prop- 
agate experimental errors, has been recently proposed [3]. 
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other words, the benchmarking role of the proton (or deuteron) cross-section in the 
experimental data is played in the global fits by a known set of free proton PDFs. 

On the other hand, the abundant neutrino DIS data available for nuclear targets, 
both iron and lead, could provide strong constrains to the global fits of nuclear PDFs. 
The absence of a corresponding proton benchmarking data in this case indicates, how- 
ever, that the price to pay is a modification of the usual procedure to check the com- 
patibility among the different data sets within a DGLAP analysis. Several attempts 
to study nuclear effects in neutrino DIS exist [101 [HI [121 [131 [H] • If has been recently 
pointed out [H], that the NuTeV neutrino- Iron DIS data [15] displays a behaviour 
that hints to a possible disagreement between the nuclear modifications Rf{x,Q'^) in 
neutrino DIS, and those extracted from charged lepton DIS and DY process. If true, 
this would seriously doubt the validity of QCD factorization for bound nucleons. This 
is not only a problem for the nuclear PDF analyses, but would also make these data not 
suitable to be included in free proton global fits, and relevant constrains for different 
flavors would be lost. 

In this work we address this issue by studying at the same time the compatibility of 
neutrino DIS data on iron and lead targets from three experiments within the approach 
of factorizable universal PDFs. These PDFs are taken from the latest available sets 
CTEQ6.6 [T] for protons and EPS09 [5] for nuclei, in which these data have not been 
included. We use, in particular, data from three collaborations NuTeV [15], CDHSW 
[I6] , and CHORUS [17] . The error analyses performed for PDFs allow us to propagate 
the uncertainties from other sets of data and to check the compatibility both between 
data and theory and among different sets of experimental data. In order to present the 
data in a more transparent manner, we provide as benchmark cross-sections for the 
nuclear case the ones computed theoretically using the known set of free proton PDFs. 
Our analysis finds a good description of neutrino DIS with nPDFs fitted to other data 
sets and do not find support for a breaking of the universality of the nuclear PDFs. We 
also comment on the compatibility between different data sets, in particular concerning 
the systematics in energy of the x-dependence in NuTeV data. 



2 Framework 



2.1 Computation Framework 



A general expression for the differental DIS cross-sections can be written down in terms 



of three nuclear structure functions _p/'^'''^ as 



2^i'A,VA 



dxdy 



(2) 



2 TjiuAjuA 

xy Fi 



i-y 



xyM\ 



^uA.uA 



± ( 1 



:ivA,VA 



2 



where y = Q^/{2xME^), Gp = 1.16637 x lO^GeV^ is the Fermi constant, M = 
0.938 GeV denotes the nucleon mass, M\y = 80.398 GeV the mass of the W boson, and 
El, to the neutrino energy. The " + " sign in front of the last term refers to the neutrino 
induced process, whereas " - " sign corresponds to the antineutrino scattering. 

Each of the nuclear structure functions are composed of those of the bound protons 
and neutrons 

7 A — 7 

jpA ^proton, A ^ ^ ^neutron, A 

We calculate the structure functions in the next-to-leading order QCD-improved parton 
model, schematically as a convolution between the perturbative Wilson coefficients u 
and the PDFs /p"*"-'^ 

^proton,A^^^^^^^proton,A^ (4) 
k 

where the k runs over all parton flavors. The neutron structure functions p^'^^^"^'^"^'^ 
computed correspondingly by assuming the isospin symmetry, i.e. 

^neutron, A _ /.proton, A /.neutron, A _ /proton, A (r\ 
Ju ~ J d ^ J d ~ Ju ■ 

The effects of the heavy quarks are taken into account by the SACOT-prescription 
[m [191 [201 [2ll [22]. This is also the scheme which was adopted by CTEQ6.6 [1] 
analysis. This set will be used here to compute the benchmark cross-sections for the 
neutrino-nucleon cross-sections. 

The nuclear modiflcations to the free proton PDFs JP™*""^ = yCTEQe.e taken 
from the recent EPS09-analysi^ [5], which relates the nuclear PDFs to the free proton 
ones by a multiplicative factor R^{x, Q"^) as 

/r°°'^(x, Q') = Rtix, Q') /r"^(x, Q'). (6) 



2.2 Experimental Input 

The experimental data on which we base our analysis are differential charged- current 
neutrino\antineutrino-nucleus cross-section^ d'^a'^'^''^^ /dxdy measured at several val- 
ues of neutrino beam energies E^, Bjorken-x, and virtuality Q^. The NuTeV collabo- 
ration [13] provides 2618 and CDHSW 1533 data points from Iron targets, whereas 
the CHORUS collaboration [17] gives 1214 data points from Lead target measurements. 
The kinematical reach of these data are illustrated by the polygons of Figure [1] The 
CTEQ analysis which is used here as reference for the proton PDFs imposes cuts in 
experimental data which will be also taken into account here. For consistency with 
this analysis and to avoid getting contaminated from higher-twist effects we restrict 

^The EPS09-analysis was performed in the zero-mass scheme. However, in the presently studied 
region of a;, we have checked that the difference between general-mass and zero-mass schemes in 
R^{x,Q'^) are small. 

■^For a review of neutrino experiments, the reader can consult e.g. Ref. |23j . 
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Figure 1: The kinematical reach of the NuTeV, CDHSW, and CHORUS neutrino data. The 
dotted blacli hues are the Icinematical borderlines of Eq. ([7]) . 

the virtuality and the final state invariant mass = + Q'^{1 — x)/x by the 
kinematical boundaries 

gL>5GeV^ Wl,>12GeV'. (7) 

These cuts are indicated by black dottet lines in Figure [T] After applying these restri- 
tions, 2041 NuTeV data points, 700 CDHSW, and 768 CHORUS data points remain. 

2.3 Target Mass and Radiative Corrections 

The two main types of additional corrections factors that should be applied before a 
meaningful comparison to the experimental data are the target mass correction, and 
the one originating from the electroweak radiation. 

The target mass corrections are important at large x and low Q^. In this work, we 
apply the collinear factorization prescription from Ref. [23], which basically amounts 
to changing the Bjorken-x variable to the Nachtmann variable, defined as ^ = 2x/(l + 
^/l + 4x^M2/Q^). To be precise, we make the substitution 

for the subprocesses involving light partons. 

In contrast to the target mass corrections, the radiative corrections are rather a 
small-x issue and more spread out in Q^. To account for the radiative corrections, we 
use the parton flavor dependent factors A^^*^'^*'^° from Ref. [25] correcting the leading 
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order terms, i.e. we compute 



k 

Both corrections modify the Q^-dependence of the calculated cross-sections. In the 
next sections, we will study the Q^- dependence of the data vs. theory ratios and and by 
examining them, one can judge whether obvious additional ingredients outside pQCD 
remain — only those regions of phase space where such behaviour is not evident, can 
be reliably used to make conclusions about the nuclear PDFs. 

2.4 Probed PDF Components 

To get a rough understanding of the behaviour of the neutrino cross-sections and their 
sensitivity to the different PDF-components, it is sufficient to consider the leading 
order expressions for the structure functions. For simplicity, we suppress the Cabibbo- 
Kobayashi-Maskawa matrix elements and assume that one is at large enough to 
forget about the heavy quark mass effects. Then, for the neutrino scattering 

F^^ = 2xF^^ = 2x{d^ + s^ + b^ + u'^ + c'^) (10) 

F^^ = 2{d^ + s^ + b^-u^-c^), (11) 

where e.g. 

= ^^proton.A ^ ^ " ^ ^ncutron,A ^3) 

7 A — 7 

_ ^proton, A _|_ ^proton, A 

^ A A 

assuming the isospin symmetry. For the antineutrino scattering 

F2^^ = 2xFf^ = 2x(m^ + c^+/ + I^ + 6^) (13) 
F^^ = 2(u^ + c^-d'^ -s^-b'^y (14) 

By inserting in Eq. (j2]), one finds that the cross-sections depend on the following 
partonic combinations: 

dV^ oc (d^ + + 6^) + (1 - y)' (n^ + c^) (15) 
dV^ oc + + +(l-y)'(M^ + c^). (16) 

Because of the suppressing (1 — factor, the antineutrino cross-section carries more 
sensitivity to the sea quarks towards larger Q"^, whereas the neutrino cross-sections 
always get a significant contribution from the valence quarks. 
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2.5 The Method Of Comparison 



In what follows, we will present the comparison to the data in terms of the following 
data vs. theory ratios: 

pCTEQ6.6 ^ g'^'^ (Experimental) 

" (T-'^ (CTEQ6.6) ' ^ ^ 

where CTEQ6.6 denotes the set of PDF used in the calculation. The PDF uncertainties 
are considered by calculating the upper and the lower cross-section uncertainties Ao"^ 
by a formula 

1 2 



k 



max 
min 



{aiS^)-a{S'),a{S^)-aiS'),0} 



where cr(S'°) is the value of the cross-section computed by the CTEQ6.6 central set, 
and c(>S'^) are the ones calculated by the kth. uncertainty sets. We will compare such 
"pseudodata" to the expected nuclear effects using the nuclear modification factors 
Rt{x,Q^) from EPS09-analysis 

OCTEQ6.6XEPS09 ^ g--'^ (CTEQ6.6 X EPS09) 

a-'^ (CTEQ6.6) ' ^ ' 

which allow us to make a conclusion about the consistency of the nuclear effects be- 
tween neutrino scattering and the charged lepton scattering. We note that the EPS09 
includes a possibility for a similar uncertainty analysis as in Eq. f lTSj) . The size of the 
uncertainties of this origin will be addressed in section 13.31 



3 Results and Discsussion 
3.1 The x^- values 

The standard way of measuring the agreement between the experimental data and the 
theory, is to calculate the defined as 

E (^)^ m 

iGdata points 

where Tj is the theoretical value, and Di is the corresponding experimental value for 
the ith. data point with estimated uncertainty (Tj 0. In this work we add the statistical 
and systematic uncertainties in quadrature when computing x^- Usually, if the x^/^ 
[N being the number of data points) is not much larger than one, the theoretical 
calculations are considered as being statistically consistent with the data. 

''The experimental data used here have been taken from the database [55]. Statistical and sys- 
tematic errors have been added in quadrature. Additional overall normalization errors have been 
ignored. 
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RAD + TM CTEQ6.6 CTEQ6.6xEPS09 
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CHORUS 
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0.71 


RAD + No TM 




CTEOfi fixEPS09 
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l.zy 


1 no 
l.Oz 


CHORUS 


1.03 


0.85 


CDHSW 


0.88 


0.66 


No RAD + No TM 


CTEQ6.6 


CTEQ6.6xEPS09 


NuTcV 


1.35 


1.08 


CHORUS 


1.23 


1.09 


CDHSW 


0.96 


0.86 



Table 1: The x^/-^- values computed using CTEQ6.6 with and without nuclear modification 
from EPS09. The numbers are given for calculations including radiative and target mass 
correction (RAD + TM), with radiative but without target mass correction (RAD + No 
TM), and without radiative nor target mass corrections (No RAD + No TM). 

In Table [H we give the x^/-^- values, computed with those data points that meet 
our kinematical constraints. We have calculated the values in three different ways 
using CTEQ6.6 PDFs with and without nuclear modifications from EPS09: With the 
radiative and the target mass corrections, with the radiative but without target the 
mass correction, and finally without the radiative nor the target mass corrections. We 
make the following observations: 

• Whatever way we make the calculation, the one without nuclear corrections from 
EPS09 gives substantially larger . 

• For CHORUS and CDHSW data, the x^s get consistently smaller as the radiative 
and target mass are applied. 

• For NuTeV data, the remains practically unchanged whether radiative or 
target mass corrections are applied or not. 

Whereas the good x^/^- values indicate that the theory calculations are in good 
agreement with the data, it may seem puzzling that while x^s for CHORUS and 
CDHSW gain improvement when applying the radiative and the target mass correc- 
tions, the XNuTeV stays essentially constant. 

The fact that the Xchorus ^^"^ Xcdhsw S^^ clearly smaller after addition of the 
radiative and target mass corrections indicates, that there are a substantial changes in 
the computed cross-sections. Therefore, the unchanging XnuTbY implies that there is 
something peculiar in the NuTeV data sample. 

We note that the target mass correction, Eq. (|H]), was not applied in the EPS09- 
analysis. Repeating the analysis including such effect would slightly change the shape 
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Figure 2: Examples of Data/Theory ratios for the CHORUS neutrino data. The black 
squares with error bars correspond to the ratios 7jCTEQ6.6^ whereas the red lines are the 
^CTEQ6.6xEPS09 predictions. The relative CTEQ6.6 uncertainty is shown by blue bands. 

of the valence quark nuclear modifications at large-x. Moreover, the CTEQ6.6 uncer- 
tainties also grow towards large x. Both effects, if taken into account, would modify 
the evaluation of the data in this region of phase space. Therefore, the changes 
in when applying the target mass corrections for the CTEQ6.6xEPS09 column of 
Table [1] should be taken only as indicative. Note, however, how the gets clearly 
worse when adding the target mass correction if the nuclear PDFs are not used. 
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Figure 3: The Q^-averaged CHORUS neutrino data. The black squares with errorbars show 
the -Ravp^S^^ ratios, while the red hnes correspond to -Raw't^S-o' 
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Figure 4: As Figure [3l but for CHORUS antineutrino data. 

3.2 Data vs. Theory Ratios 

In this section we take a closer look to the nuclear effects in neutrino-nucleus DIS by 
plotting results for ratios i?CTEQ6.6 ^^^^ ^cteq6.6xepso9 ^^^^^ ^j^j^^ ^j^g Figure M 

shows two representative examples of the energy and Q^-dependence of the data vs. 
theory ratios. Evidently, as one cannot observe a large systematic Q^-dependence in 
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Figure 5: As Figure [3l but for CDHSW neutrino data. 
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Figure 6: As Figure El but for CDHSW antineutrino data. 



these panels, it is reasonable to look a Q^-averaged versions of these figures in order to 
better see the big picture — the x-shape of the nuclear modifications. 
We will form the Q^-averages of the panels by the prescription 



pCTEQ6.6 ^ 
-""Average — I / j 



( ^ j^^CTEQe.i 



iGfixcda; / ViSfixedx 



(21) 



\iGfixod X 

where the sum runs over all data points in same x- and energy-bin, but different Q^. We 
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Figure 7: As Figure El but for NuTeV neutrino data. 



note that the way we compute the average and, especially its uncertainty, is not meant 
to be statistically "the correct way" , but it is rather ment to give an idea about the size 
of the uncertainties. We stress, however, that the statistical analysis (the calculation of 
the x^-values in the previous sections) was performed with the original (non-averaged) 
data. The purely theoretical calculations are simple averages, for example 



pCTEQ6.6xEPS09 _ J_ pCTEQ6.6xEPS09 /r,9\ 

-"-Average — / ^ -'^i ■ 

jg fixed X 



The results of this exercise are presented in Figures |3]|H1 We wish to draw the reader's 
attention to the following points: 

• The CTEQ6.6xEPS09-prediction and the CHORUS and CDHSW neutrino data 
are in a perfect agreement. The data do not display an evident dependence of 
the neutrino energy. This is also true for the CDHSW antineutrino data. 
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Figure 8: As Figure El but for NuTeV antineutrino data. 



• The CHORUS antineutrino data seems to somewhat overshoot the CTEQ6.6xEPS09- 
prediction around x ~ 0.2 staying, however, clearly in the uncertainty limits. 

• Whereas some panels corresponding to the NuTeV neurino data are very well in 
line with the CTEQ6.6xEPS09 prediction, there is an evident variation from a 
beam energy to another. For example Ei^ = 170 GeV and E^y = 190 GeV are in 
the brink of disagreement. 

• Similar evident energy-dependece is not observable in the NuTeV antineutrino 
panels. Within the uncertainties the data are well reproduced. 

The clear energy- dependent panel-to-another variation in the NuTeV neutrino data/theory 
ratios is evidently something that cannot be compensated by only re-fitting the PDFs 
— the CTEQ6.6xEPS09 predictions vary very slowly as a function of the incident neu- 
trino energy. On the contrary, using this data set as a source of PDFs in a fit would 
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inevitably result in a compromise between data in different energy bins displaying mu- 
tual tension. No signs of this kind of evident controversy are visible in the NuTeV 
antineutrino data nor in the CHORUS or CDHSW data. 

This observation also explains why the XnuTcV almost inert to the target mass and 
radiative corrections. Indeed, by simply ignoring the Ei, = 35, 45, 130, 150, 170, 245 GeV 
panels when computing XnuTgY' would cause the XnuTov/^ ^'^ reduce from 1.10 to 0.94 
when applying the corrections. 

3.3 Other Sets of nPDFs 
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Figure 9: Two panels of Q^-averaged CHORUS and NuTeV data with predictions from three 
sets of nPDFs, EPS09 (black line with blue errorband), HKN07 (red dashed line), and nDS 
(blue dotted line). 



So far we have considered the nuclear effects in PDFs utilizing solely the latest 
available analysis, EPS09, observing a good average agreement with the experimental 
data. To address whether the same conclusion would have been drawn using different 
sets of nPDFs we show, in Figure [9l examples of the Q^-averaged ratios from the two 
other available NLO sets nDS and HKN07, and compare them to the uncertainty band 
of EPS09. 

The upper panels illustrate the situation with the neutrino cross-sections. At large 
X, the dominant contribution to the cross-sections comes from the valence quarks and 
— due to good amount of valence quark constrains from charged lepton DIS there — 
the different sets of nPDFs give predictions relatively close to each other. Only the 
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nDS prediction for Lead display a clear deviation from others. Towards small x, the 
predictions from different sets of nPDFs are very close to each other. 

In the case of antineutrino cross-sections, illustrated in the lower panels, there is an 
enhanced contribution from the sea quarks, for which the nuclear modifications at large 
X are not very well constrained. This uncertainty is reflected as clearly larger scatter 
between the curves. At small x, all sets of nPDFs are again in a good agreement. 

In summary, it is mainly in the large-x sector where the differences between sets 
of nPDFs become apparent: While the nDS description admits some improvements 
especially regarding the A-dependence in this region, the HKN07, on the other hand 
clearly misses the structure of antineutrino cross-sections at large-x. This is due to 
the radically different behavior of antiquarks and gluons when compared with the two 
other sets in the corresponding relevant region of phase space. These findings point 
also to the adequacy of neutrino data in providing further constrains for nPDFs. 

4 Conclusion 

In this article we have presented a study of the nuclear effects in the (anti)neutrino- 
nucleus deeply inelastic scattering. As a conclusion, we find that from the analysed 
data sets only one — the NuTeV neutrino data — shows signs of not being in a 
complete agreement with the present-day PDFs. This discrepancy, we argue, is however 
something that cannot be completely cured by simply re-fitting the PDFs. Rather, it is 
due to an unexplained, neutrino energy-dependent fluctuations in the data — in order 
to make this point more explicit, we have presented Q^-averaged data where the energy 
systematics of the x-dependence is more transparent. Therefore we assert that this 
data set cannot be taken as a discriminating factor when making conclusions about 
universality of the nuclear PDFs. The forthcoming data from NOMAD experiment 
will eventually facilitate a further clarification about the conflict between the different 
experiments. 

Apart from this finding, the present analysis do not give any reason to belive that 
the nuclear effects in neutrino-nucleus DIS would be different from those extracted 
from charged lepton DIS and DY production of dileptons — the factorization seems to 
work well. 

Having said that, there is a lot of room for improvement and future work. For 
example, the EPS09-analysis leaves no freedom to the strange quark nuclear modifica- 
tion, but it is something that is fixed by hand to other sea quarks. As the antineutrino 
cross-sections carry a non-neglibible contribution from the strange quarks, including 
these data to a global fit of nPDFs would offer a chance to free this assumption made 
in EPS09. This could improve the slight undershooting of the CHORUS antineutrino 
data. Also, whereas the Iron is almost an isoscalar nucleus {Z = 26, = 30), the Lead 
is clearly non-isocalar [Z = 82, = 126). This could offer a possibility to separately 
constrain the nuclear modifications for the up and down quarks which has not been 
possible before. 
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